Roughening transition in a model for dimer adsorption and desorption 
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A solid-on-solid growth model for dimer adsorption and desorption is introduced and studied nu- 
merically. The special property of the model is that dimers can only desorb at the edges of terraces. 
It is shown that the model exhibits a roughening transition from a smooth to a rough phase. In 
both phases the interface remains pinned to the bottom layer and does not propagate. Close to the 
transition certain critical properties are related to those of a unidirectionally coupled hierarchy of 
parity-conserving branching-annihilating random walks. 



PACS numbers: 64.35.Rh, 05.70.Ln, 82.20.Wt. 

The study of crystal growth and transitions between 
different morphologies of moving interfaces is a field that 
continues to attract great interest |ffl]. Various models 
have been developed in order to describe the essential 
features of roughening transitions . In many cases the 
critical properties close to the transition are universal, 
i.e., they do not depend on the microscopic details of 
the model. For example, the majority of roughening 
transitions in d > 2 spatial dimensions belongs to the 
Kardar-Parisi-Zhang universality class ||. Therefore it 
is an important theoretical task to categorize the possible 
universality classes for roughening transitions in a given 
geometry. 

It is well known that one-dimensional interface mod- 
els with short-range interactions at thermal equilibrium 
do not exhibit roughening transitions. However, under 
nonequilibrium conditions such a transition may occur 
even in one spatial dimension. Known examples in- 
clude polynuclear growth (PNG) models [Q, solid-on- 
solid models with evaporation at the edges of terraces || , 
and certain models for fungal growth ||. The key fea- 
ture of all these growth processes is a close relationship to 
directed percolation (DP) [Rj^j. More precisely, the DP 
process emerges at a particular reference height of the in- 
terface. In the models of Ref. || this reference height is 
the spontaneously selected bottom layer of the interface, 
whereas in polynuclear and fungal models, which evolve 
by parallel dynamics, the reference level propagates at 
maximal velocity. The sites where the interface touches 
the reference height correspond to the active sites of DP. 
Therefore, in the active phase of DP the interface fluc- 
tuates close to the reference level so that the interface is 
smooth. On the other hand, in the inactive phase of DP, 
the interface detaches from the reference level and evolves 
into a rough state. More recently it has been shown that 
the critical behavior at the first few layers next to the 
reference height is equivalent to that of unidirectionally 
coupled DP processes || . 

DP itself is the generic universality class for phase tran- 
sitions into absorbing states and covers a wide range of 
models. Only a few exceptions from DP are known. 



One of them is the so-called parity- conserving (PC) 
class |]lo|-|l5| which is represented most prominently by 
branching-annihilating random walks with even number 
of offspring (BAWE) 0. The PC class also includes 
models with two symmetric absorbing states p^ , [l4] . |l5| . 
where the kinks between different absorbing domains 
may be interpreted as walkers whose number is preserved 
modulo 2. In other models Jl(],[l3| the two symmetric ab- 
sorbing states emerge as checkerboard-like configurations 
of particles at even or odd sites, respectively. 

Recently Park and Kahng jl6| posed the question 
whether it is possible to replace the underlying DP mech- 
anism in the interfacial growth models of Refs. by 
a PC mechanism. To this end they introduced a model 
which involves two symmetric particle species and ex- 
hibits a roughening transition in 1+1 dimensions. Since 
the dynamical rules at the bottom layer mimic a con- 
tact process with two symmetric absorbing states [ f[5| , 
the authors expected the roughening transition to be re- 
lated to the PC class, in the same way as the transition 
in the monomer models of Ref. j5j is related to DP. But 
surprisingly the numerical estimates for the critical ex- 
ponents deviated significantly from the PC values. Park 
and Kahng argued that the unexpected behavior may be 
related to the fact that dynamical processes at lower lev- 
els are strongly suppressed by particles at higher levels. 
In particular, kinks between different particle species may 
become frozen when they are covered by another layer of 
atoms. 

In this Letter we introduce a model for dimer adsorp- 
tion and desorption where this suppression effect is much 
less important. As dimers consist of two atoms, the num- 
ber of particles at each height level is preserved modulo 2. 
The dynamical rules are defined in a way that they mimic 
a BAWE at the bottom layer of the interface. It turns 
out that some of the critical properties of the roughening 
transition in this model are indeed characterized by PC 
exponents. 

a. Definition of the model. The model is defined on 
a d-dimensional square lattice with L d sites and periodic 
boundary conditions. The interface height at site i is 
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FIG. 1. Dynamical rules: Dimers are adsorbed with prob- 
ability p and desorbed at the edges of terraces with probability 
1-p. 



represented by an integer height variable hi which obeys 
the restricted solid-on-solid (RSOS) condition 
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The model evolves by random sequential updates accord- 
ing to the following dynamical rules (see Fig. |l|). In each 
attempted update two adjacent sites i\ and i% are se- 
lected at random. If the heights hi 1 and hi 2 are equal, 
one of the following moves is carried out. Either dimers 
are adsorbed with probability p 

hi ± —yfiij^ + 1, h i2 — ► h i2 + 1 , (2) 

or dimers desorb with probability 1 — p 
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where min^/^^ hk denotes the minimum height of sites 
i\ and ii and their nearest neighbors. An update will be 
rejected if it leads to a violation of the restriction (|l|). 

The above model is defined in arbitrary spatial dimen- 
sions d and is translationally invariant in space, time as 
well as in height direction. Moreover, it may be eas- 
ily generalized to rt-particle objects (monomers, dimers, 
trimers, etc.). In particular, for n = 1 it reduces to 
the previously studied monomer growth model of Ref. |^] 
which is related to DP. We use random sequential up- 
dates in order to demonstrate that in contrast to PNG 
models QJ^] , the phase transition in our model does not 
require synchronous updates. We have also studied an- 
other variant of the dimer model with parallel updates, 
as will be described elsewhere |l7| . This cellular automa- 
ton model resembles a layer-by-layer crystal growth pro- 
cess and exhibits the same universal properties as the 
random-sequential variant. 

b. Phenomenology of the roughening transition. If the 
adsorption rate p is very small, only a few dimers are 
adsorbed, staying for a short time, before they quickly 
evaporate again. Therefore, the interface is smooth and 
anchored to the bottom layer (the spontaneously selected 
line of minimal height). As p increases, more and more 



dimers cover the surface and large islands with several 
layers stacked on top of each other are formed. When p 
exceeds a certain critical value p c , the mean size of the 
islands diverges and the interface heights become uncor- 
rected over long distances, i.e., the interface evolves into 
a rough state. 

In analogy to the monomer models of Ref. @ we may 
expect the interface to detach from the bottom layer in 
the rough phase, resulting in a finite propagation veloc- 
ity. However, it turns out that in the present case the 
interface remains pinned to the initial height and does 
not propagate at constant velocity. This is due to the 
fact that a stochastic deposition process cannot create 
a dense packing of dimers. The emerging configurations 
are rather characterized by a certain density of defects 
(solitary sites at the bottom layer) where dimers cannot 
be adsorbed. Because of the RSOS condition (|l|) these 
defects act as 'pinning centers' which prevent the inter- 
face from growing. The pinning centers cannot disappear 
spontaneously, rather they can only diffuse by interface 
fluctuations and recombine in pairs so that their number 
is expected to decrease extremely slowly. Therefore, the 
interface of an infinite system in the rough phase does 
not propagate at constant velocity. Instead the average 
height increases only logarithmically with time, as will 
be shown below. 

In order to illustrate the PC mechanism in this model, 
let us interpret the sites at the bottom layer hi = as 
A-particles. Adsorption and desorption processes corre- 
spond to certain effective reactions of the A-particles. For 
example, the adsorption of dimers at the bottom layer 
corresponds to a pair-annihilation process 2A — > at 
rate p. On the other hand, when dimers evaporate, two 
A-particles are created. However, since dimers can only 
evaporate at the edges of terraces, the presence of an- 
other neighboring A-particle is required, giving rise to 
an effective reaction A — > 3A at rate 1 — p. These two 
processes compete one another and realize a BAWE. This 
mapping is, of course, not exact since the A-particles are 
also coupled to the dynamical processes at higher levels 
of the interface. Nevertheless it seems that this feedback 
does not affect the critical behavior of the bottom layer. 

c. Critical properties of the interface width. The or- 
der parameter that describes a roughening transition is 
the interface width 
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In order to investigate its critical properties, we perform 
Monte-Carlo simulations in d — 1 dimensions on a lattice 
with L — 4 096 sites up to almost 10 5 time steps. Start- 
ing from an initially flat interface hi = we find the 
following results. Below the critical threshold p < p c the 
width saturates at some finite value which means that 
the heights hi are correlated over long distances. At the 
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FIG. 2. Critical properties of the interface. In parts (a) 
and (b) the (squared) width and the average height are shown 
as functions of time near and far away from criticality, respec- 
tively. Parts (c) and (d) show data collapses for finite-size 
simulations (see text). 



roughening transition p — p c ~ 0.317(1) we observe that 
the squared width and the average height grow logarith- 
mically with time as (see Fig. 0a) 



W 2 {t) ~h(t) ~ log* 



(5) 



In the rough phase p > p c this temporal behavior crosses 
over to W(t) ~ log* (see Fig. |]b). 

Assuming that the roughening transition is governed 
by an underlying PC transition, it is near at hand to 
conjecture the finite-size scaling form 



W 2 (L,t) - log r a F(t/L z ) 



(6) 



where F is an universal scaling function and z = v\\/v±_ = 
1.76(5) is the dynamical scaling exponent of the PC class. 
The exponent a is unknown and plays the role of a rough- 
ening exponent. In order to verify this scaling form, 
we perform finite-size simulations at criticality and plot 
exp(W 2 )/t a for 10 2 < t < 10 5 against t/L z . Using the 
estimate a = 0.172(10) we obtain a fairly accurate data 
collapse (see Fig. ||c), supporting that the dynamical ex- 
ponent in this model is the same as in PC transitions. 

d. Critical properties of the first few layers. We now 
turn to the critical behavior at the first few layers of the 
interface. Let us denote by 
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FIG. 3. (a) Densities n k at the first three levels as a func- 
tion of time at and close to criticality. (b) Particle densities 
in unidirectionally coupled BAWE's at criticality. 

the density of sites i with hi < k. In analogy to Ref. 0] 
we expect n k to obey the finite-size scaling form 



n k (L,t)~t~ s »G(t/L z ) , 



(8) 



where G is a universal scaling function. So = S = /3/Vii 
denotes the usual cluster survival exponent of the PC 
class while the other exponents <5i,<52, ••• should be in- 
dependent and smaller than <5o- Starting from a flat in- 
terface at h = in a large system with 4 096 sites, we 
measure no, n\ and ri2 at criticality (see Fig. |3^,) which 
should decay as nk(t) ~ t~ Sk . Averaging the slopes over 
the last decade in time we obtain the estimates 

<5 = 0.280(10) , <5i = 0.200(15) , S 2 = 0.120(15) . (9) 

The estimates for <5 and z can also be verified by a finite- 
size data collapse (see Fig. ||d) . The numerical value of <5 
is in fair agreement with the PC exponent 5 = 0.285(5), 
confirming that the dynamical processes at the bottom 
layer belong to the PC universality class. 

e. Unidirectionally coupled BAWE's. In order to ex- 
plain the exponents 5i and 62 we propose that the crit- 
ical behavior at the first few layers in the dimer model 
should be equivalent to that of a hierarchy of unidirec- 
tionally coupled PC processes, just in the same way as 
the monomer model of Ref. E] is related to unidirection- 
ally coupled DP processes £]. To this end we extend 
the particle interpretation: assuming that h = is the 
bottom layer of the interface, let us interpret sites i with 
hi < 0,1,2,... as particles A, B, C, . . ., respectively. By 
definition, particles of different species are allowed to oc- 
cupy the same site simultaneously. As shown before, the 
temporal evolution of the A-particles resembles a BAWE. 
Similarly, the B-particles perform an effective BAWE on 
top of inactive islands of the ^4-system. Clearly, the tem- 
poral evolution of the i3-particles depends strongly on the 
dynamics of the A-particles. For example, A-particles in- 
stantaneously create i?-particles at the same site, which 
may be interpreted as an effective reaction A — > A + B 
at infinite rate. As this reaction does not modify the 
configuration of the A-particles, it couples the two sub- 
systems only in one direction without feedback. On the 
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FIG. 4. Unidirectionally coupled BAWE's. 



other hand, the RSOS condition (Q) introduces an effec- 
tive feedback so that the A-particles are not completely 
decoupled from the £>-particles. However, we will assume 
that in contrast to Ref. |ll| this feedback is irrelevant for 
the critical behavior and can be neglected. Similarly, the 
C-particles are coupled to the B-particles by the effective 
reaction B — > B + C. Therefore the dimer model resem- 
bles a linear sequence of BAWE's which are effectively 
coupled without feedback in one direction, as sketched in 
Fig. |. 

To verify this conjecture, we study three unidirection- 
ally coupled BAWE's by Monte Carlo simulations (de- 
tails will be published in Jl7[] ) . We implement the cou- 
pling between the systems in a way that A-particles in- 
stantaneously create B-particles at the same site in the 
-B-system, provided that this site is empty. Similarly, 
B-particles create C-particles at the same site. We sim- 
ulate a hierarchy of three coupled BAWE's on a lattice 
of 2 500 sites up to 25 000 time steps and measure the 
densities of A, B, and C-particles as functions of time av- 
eraged over 1000 independent realizations (see Fig. ||b). 
It turns out that the particle densities tla , tib , nc de- 
cay in the same way as no,ni,ri2 in the dimer model. 
Averaging the slopes over the last decade in time we ob- 
tain the exponents 6 A = 0.280(5), S B = 0.190(10), and 
5c = 0.120(15), which are in fair agreement with the 
corresponding exponents of the dimer model in Eq. <^). 
This supports the hypothesis that the dimer model and 
coupled BAWE's exhibit the same type of critical behav- 
ior. 

It should be noted that the densities at higher levels do 
not scale perfectly. Rather the curves for k > 2 in Fig. || 
appear to be slightly bent, indicating a deviation from 
perfect scaling. This curvature is neither related to errors 
in the estimation of p c nor to finite-size effects. Similar 
deviations were also observed in the case of DP-related 
growth models |9) and are presumably caused by certain 
infrared divergent diagrams of the underlying field the- 
ory. 

To summarize, we have introduced an interface model 
for dimer adsorption and desorption and have presented 
numerical evidence that the roughening transition in this 
model may be associated with the universality class of 
unidirectionally coupled PC processes. 
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